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0\ ' Abstract 

,-J ' We point out that the cosmological moduli problem is not necessarily resolved even if the 

(T) ■ modulus mass is heavier than 0(10) TeV, contrary to the common wisdom. The point is that, 

in many scenarios where the lightest moduli fields are stabilized by supersymmetry breaking 

K> | effects, those moduli fields tend to mainly decay into almost massless axions, whose abundance 

H ■ . 

is tightly constrained by the recent Planck results. We study the moduli-induced axion problem 

in concrete examples, and discuss possible solutions. The problem and its solutions are widely 
applicable to decays of heavy scalar fields which dominate the energy density of the Universe, 
for instance, the reheating of the inflaton. 



I. INTRODUCTION 

The cosmological moduli problem [l|, |2| is one of the most important issues in string 
theory and cosmology. The purpose of this paper is to point out that the cosmological 
moduli problem is not necessarily resolved even if the modulus mass is heavier than 
O(10) TeV, contrary to the common wisdom. 

In superstring theories |3j, the so-called moduli fields appear at low energies through 
supersymmetric (SUSY) compactifications e.g., on a Calabi-Yau (CY) space |4[ as Kaluza- 
Klein zero-modes of 10-dimensional metric and p-form fields. These moduli and their 
axionic superpartners are massless at the perturbative level because of the shift symmetry; 

-'moduli ' -'moduli ' 10L-, \i-) 

which is regarded as a remnant of higher dimensional gauge symmetry, with a being 
a real transformation parameter. In order to have a sensible low-energy theory, those 
moduli fields must be stabilized. It is well known that many moduli fields are fixed 
simultaneously by the closed string flux backgrounds in extra dimensions, i.e. flux com- 
pactifications p, |6[, and most of the remaining moduli not fixed by the fluxes can be 
stabilized by instantons/gaugino condensations a la KKLT |7|. 

During inflation, some of those moduli fields, especially relatively light ones, are likely 
deviated from the low-energy minima. They will start coherent oscillations with a large 
amplitude at some time after inflation, and soon dominate the energy density of the 
Universe. If its mass is of order the weak scale, it typically decays during the big bang 
nucleosynthesis (BBN), thus spoiling the overall agreement between BBN theory and light 
element observations. If the modulus mass is much lighter, the situation becomes worse; 
the modulus abundance can easily exceed the dark matter abundance or its decay may 
produce too much X-rays or gamma-rays |8j. This is the notorious cosmological moduli 
problem l| . Among many solutions proposed so far, the simplest one is to assume a heavy 
modulus mass; the BBN bound is relaxed significantly if the modulus mass is heavier than 



several tens TeV. 1 

The heavy moduli scenario does not necessarily lead to a successful cosmology. It 



was pointed out in Refs. [HH12J that gravitinos are generically produced by the modulus 
decay if kinematically allowed, and that gravitinos thus produced affect the light element 
abundance and/or produce too many lightest supersymmetric particles (LSPs), even if 
the modulus decays before BBN. This is known as the moduli-induced gravitino problem. 
There are several ways to avoid the problem. For instance, no gravitinos are produced if 
the decay is kinematically forbidden. This requires a heavy gravitino mass comparable 
to or heavier than the modulus. Alternatively, even if many gravitinos are produced by 
the modulus deca y, t he cosmological bound can be relaxed if the gravitino is heavier than 
several tens TeV 13J and if the R-parity is broken by a small amount. 

There exists yet another serious cosmological obstacle. To see this, first let us note that 
both real and imaginary components of the moduli fields acquire the same mass, if they are 
stabilized by a large SUSY mass. This is the case if the moduli are stabilized by the fluxes 
or by the KKLT mechanism. On the other hand, some of the moduli may be stabilized 
by SUSY breaking effects such that their axionic fields remain extremely light due to 
the shift symmetry ([I]). Indeed, in many string models, there are often such ultralight 



axions 



14 



2l|. Those axions remain massless unless the shift symmetry is broken by an 
appropriate non-perturbative effects generated in the low energy. Furthermore, their real 
component partners tend to be lighter than those stabilized a la KKLT, because they are 
stabilized through the SUSY-breaking effect and their masses are comparable to or lighter 
than the gravitino mass. Therefore, it is crucial to study the cosmological impact of such 
light moduli stabilized by SUSY breaking effects, as the lightest moduli fields usually play 
the most important role in cosmology. 

In this paper we will show that such modulus generally decays into a pair of ax- 
ions, contributing to dark radiation whose abundance is tightly constrained by the recent 



Planck data 



22j. We call this problem the moduli-induced axion problem. As we shall see 



Sec Rcfs. M for other solutions. 



shortly, this constrains a large portion of the parameter space, and most importantly, the 
problem persists even for a modulus mass heavier than several tens TeV. The presence 
of such light moduli and ultralight axions may be a natural outcome of the string theo- 
ries, although it certainly depends on the details of the model such as the properties of 



compact geometry and brane configurations, 
by the string theoretic QCD axion 
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ndeed, if the strong CP problem is solved 



20] , it implies that there is at least one such 



modulus. Therefore, we believe that the moduli- induced axion problem is universal, and 
its solutions will provide us with important information on the high energy physics. 

L astly let us mention the related works in the past. It was recently pointed out in 
Ref. 23j that the modulus decays into a pair of its axions in a context of the (moder- 
ately) LARGE volume scenario (LVS) 24j, and the produced axions will behave as extra 



radiation since the axions are effective 
been extensively studied in Refs. [25 



y 



massless. Furthermore, the decay process has 



26j, focusing on a possibility that the produced 



axions explains the excess of dark radiation hinted by the observation at that time. In 
the context of the SUSY QCD axion, it is well known that the QCD saxion tends to decay 
mainly into a pair of QCD axions 27] . The abundance of relativistic axions produced by 
the saxion decays was studied in Refs. 28H35J . In particular, several ways to suppress the 
branching fraction of the saxion decaying into axions were discussed in Ref. 35]. 

The rest of this paper is organized as follows. In Sec. [Ill we evaluate the modulus decay 
rate into axions and the Standard Model (SM) particles in a general setup and show the 
robustness of the moduli- induced axion problem. Possible solutions to the problem will 
also be mentioned. In Sec. IHH we study concrete supergravity setups motivated by the 
string theory in order to illustrate the moduli-induced axion problem. Sec. IIVI is devoted 
to conclusions and discussion. 



II. MODULI-INDUCED AXION PROBLEM 

Let us start with a simple low-energy effective theory containing one light modulus 
stabilized by SUSY breaking effects. We consider the following Kahler potential, 

K(T,T*) = K(T + T^), (2) 

which respects the shift symmetry (CQ). Here and in what follows, we call the real com- 
ponent of T the modulus, whereas the axion refers to its imaginary component. For the 
moment we assume that the modulus is stabilized by the Kahler potential and it does not 
appear in the superpotential so that axion remains massless. Later we consider a case 
where the axion has a small but non-zero mass. The relevant terms in the Lagrangian 
are summarized in Appendix |AJ. In the following we adopt the Planck unit in which 
Mp ~ 2.4 x 10 18 GeV is set to be unity. We also denote the lowest component of a 
superfield by the same letter. 

Let us define 

m (m i I t + ia 

T - <T)s 7s?wr- (3) 

where r and a are (canonically normalized) real and imaginary components of T, re- 
spectively. Here and in what follows, the subscript T denotes the partial derivative with 
respect to T. The partial decay rate into a pair of axions is given by 



r« = T(r^aa)= ^#^4, (4) 



64tt K? 



where m T is the mass of r. The modulus also couples to the axino, a, the fermionic 

partner of the axion. The partial decay rate of the modulus into the axino pair is given 

by 

N 1 K 2 TTT 2 / 4m|\ 3/2 

where m a denotes the axino mass. (See Appendix |A"1 ) The rate © is suppressed by a 
factor of ~ (ma/m T ) 2 with respect to (J3J). 



2 See also \3 



In general, the modulus r decays also into matter fields. Suppose that the modulus 
has the following coupling, 

K D Z U \H U \ 2 + Z d \H d \ 2 + g(T + T^) (H u H d + h.c.) : (6) 

where H u and H d are up- and down-type Higgs superfields, respectively, Z u ^ is a Kahler 
metric for H u ^ that depends on the moduli, and g is some function of T + TK The partial 
decay rate of r into Higgs bosons is 3 

T(r^HH) ~ -L 9 l m% (7) 

87r K TT Z u Z d 

where we have neglected the mass of the Higgs bosons. The decay rate into higgsinos is 

given by 

r( T -> HH) = i- ^f m r , (8) 

Sit K TT Z u Z d 

where 

(9) 
Here, 7713/2 denotes the gravitino mass and _F T = —e K ^ 2 K T: '(DjW)* is the modulus in- 
terna. 

Similarly, if T contributes to the SM gauge kinetic function / vis , the modulus can decay 
into gauge bosons with the rate, 

r(r -> u> = i^BST^ do) 

1287T {ReJ vis ) z K TT 

where N g represents the number of gauge bosons, and it is given by N g = 8, 3, 1 
for SU(3),SU(2) and U(l), respectively. The decay into gauge bosons is sizable if 
|^r/vis|/(Re/vi S ) is of order unity. The decay rate into gauginos (A) is given by 



iT*\ 



N g I (dr/vta) (if + FT) + (d 2 T f vis )F T - 2(d T f vis )m x 



m T 



r ( T ^ AA ) :: T^ m 7 x 2 T^ (11) 

1287T (Re/ v i S )^ Att 



3 If the decay into heavy Higgs is kinematically forbidden and the modulus decays only into hh, ZZ 
and WW, the rate should be multiplied with (sin 2 2/3) /2, where tan/3 = (H° a ) / (H%). 



where m\ denotes the gaugino mass. It also depends on the modulus F-term and it is (at 
most) comparable to that into gauge bosons. In particular, for a generic Kahler potential, 



the partial decay rate into gauginos is not suppressed by the gaugino mass 371 ]. 

For a generic Kahler potential the decay rate into axions, Y{r — > aa), is comparable 
to what is expected based on the dimensional argument. Therefore there is no a priori 
reason to expect that the decay into axions is negligibly small with respect to the other 
decay processes, and so, the branching fraction is generically sizable, B a = Br(r — > aa) = 
0(0.1). 

In order to get the feeling that the branching fraction of the axion production tends to 
be large, let us consider a simple example before continuing further. In the next section 
we will study a few examples based on more realistic moduli stabilization. Consider the 
following Kahler potential of the no-scale form, 



K = -3 log 



Tt - 1 {\H U \ 2 + \H d \ 2 + (zH u H d + h.c.)} 



(12) 



= -3 log (T + Tt) + ^-^ (\H U \ 2 + \H d \ 2 + (zH u H d + h.c.)) + • • • , (13) 

where we assume that the modulus r is stabilized by higher order terms not shown here, 
and z denotes a coupling constant. The superpotential and the gauge kinetic function 
are assumed to be irrelevant for the modulus decay. The decay into the higgsino pair 
vanishes in the no-scale model as can be easily checked by using OH]). In the no-scale model 
T dominantly breaks SUSY and a becomes goldstino 'eaten' by gravitino. Since m T < 
m 3/ / 2 in the no-scale model, the decay into a pair of a (or the gravitino) is kinematically 
forbidden. Then the branching fraction of the axion production is given by 

B a = —^ . (14) 

2z 2 + l v ; 

Therefore, B a is indeed of order 0.1 for z = 0(1). 

As mentioned in the Introduction, the presence of additional relativistic degrees of 

freedom is tightly constrained by the Planck results. The constraint on the effective 



number of neutrinos, iV c fj, reads 22 ] 



N cS = 3.30toil (95%; Planck + WP + highL + BAO). (15) 



In the present scenario, AA r c fj(= N e s — 3.046) is related to B a as 

ANeS = « ( jaw y _?:_, (16) 

where g* counts the relativistic degrees of freedom at T = T rf , and T rf is the decay tem- 
perature of the modulus defined by 



1 (n 2 g 



T d = (1 - B a )i i-^j VAotaiMp, (17) 

with r total being the total decay rate of the modulus when the modulus dominates over 
the energy density of the Universe. Thus, the branching fraction is bounded above as 
B a < 0.12 ~ 0.22 for g* = 10.75 ~ 106.75. 

We show the cosmological constraints on the partial decay rates T a and Tsm i n Fig. HJ 
where Tsm denotes the decay rate into the SM particles. In the lower left shaded (green) 
region, the modulus decay temperature is lower than 6 MeV, and the 4 He abundance is too 



large to be consistent with observations [38|. 4 In the upper left shaded (pink) region, the 
axionic dark radiation is overproduced, leading to AiV c fj > 0.84. It is worth noting that 
the dark radiation constraint extends to the decay temperature much higher than 6 MeV. 
In other words, the moduli-induced axion problem is not solved by simply increasing the 
modulus mass. Note also that the constraint can be much more severe if the axion has a 
small but non-zero mass and decays into photons, electrons, etc. at late time. 

We have here neglected the modulus decay into higgsinos and gauginos because the 
rates depend on the modulus stabilization mechanism and the relation between the mod- 
ulus and higgsino/gaugino masses. Even if these decay processes are taken into consid- 
eration, the constraints shown in Fig. [T]will not be changed much. In this respect, the 
moduli-induced axion problem is robust. Note that the abundance of the LSPs produced 
by the modulus decay may exceed the observed dark matter abundance, if T^ < 1 GeV and 
R-parity is conserved. This is especially the case if the decay rate into higgsinos and/or 
gauginos is comparable to that into Higgs and/or gauge bosons. The LSP overproduction 



4 Note that our definition of T^ is higher than that defined in 38| by a factor of y/i. 
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FIG. 1: The cosmological bounds on the partial decay rates, T a and T$m, are shown. In 
the upper left shaded (pink) region, the axionic dark radiation is overproduced, leading to 
AiVeff > 0.84. In the lower left shaded (green) region, the modulus decay temperature is lower 



38J. The 



than 6 MeV, and the 4 He abundance is too large to be consistent with observations 

dotted (blue) lines are contours of the decay temperature, T^ = 4, 10, 30, 100 MeV from left to 

right. In this figure, we have not considered the LSP overproduction through the decay. 



can be avoided if the LSP annihilation cross section is relatively large or the R-parity is 
violated by a small amount. 

Lastly let us discuss possible solutions to the moduli-induced axion problem. There 
are basically two ways to solve the problem. One is to increase the decay rate into the 
SM particles by introducing additional decay channels or setting the coupling constants 
larger. For instance, in the aforementioned example, the branching fraction B a satisfies the 
Planck constraint if z > 2 (see Eq. (TT4"1)). Alternatively, we may introduce more than four 
pairs of Higgs doublet with a similar coupling to the modulus with z ~ 1. The decay rate 
into SM gauge bosons is also sizable if the gauge kinetic function is modulus dependent. 



Actually, if the SM branes are wrapping on the cycle in stringy compactifications, the 
corresponding moduli couple to the gauge bosons through gauge kinetic functions. An 
well-motivated example is the string theoretic QCD saxion, as we will see in Sec. IIIIBI 
The other solution is to introduce an approximate Z 2 symmetry on T, which suppresses 
Kttt/Ktt- The Z2 symmetry can be just a coincidence because one order of magnitude 
suppression is sufficient to satisfy the bound. However, one needs to make sure that the 
decay rate into the SM particles should not be similarly suppressed. In the next section 
we will take up two examples based on concrete moduli stabilization to illustrate the 
moduli-induced axion problem and its solutions. 

III. EXAMPLES 

In this section, we shall consider two examples, where an axion will appear at low 



energy scales: LVS [24J and KKLT-like one [lq |. In the LVS case, the axion remains 
light because large extra dimensions suppress any shift-symmetry breaking effects, while 
in the latter case there exists a specific geometry which preserves the shift symmetry. 
One will find that the production of axion dark radiation is unavoidable in the LVS once 
the lightest modulus dominates over the energy density of the Universe. On the other 
hand, in the latter case, the decay rate into axions can be suppressed by choosing an 
appropriate internal geometry, which implies a mild fine-tuning between the intersection 
number among Kahler forms on a CY space (Kahler potential) and such a geometry 
(superpotential). 
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A. LARGE volume scenario 



Let us consider an effective action of Kahler moduli on a Swiss cheese CY space in flux 



vacua 5 [241 ]: 



^moduli = -2 log (v + |Y V = (T b + T fe f ) 3 / 2 - (T s + Tlf'\ (18) 

27T 

H/ modu u = Wo + Ae- a * T °; a s = — . (19) 

Here, T& and T s are the Kahler moduli and V is the CY volume; £, H^o and A are of 0(1) 
constants. iV is a positive integer. We have assumed that the 4-cycle S supported by 



Re(T s ) is rigid divisor, i.e. h 1,0 (S) = h 2,0 (S) = 41] : A non-perturbative effect can exist 
on such a divisor. At the minimum of the scalar potential V = e K [\DW\ 2 — 3\W\ 2 ] + jD 2 , 
where DW = (dK)W + dW, one finds 

V ~ -^exp^ 3 ) » 1, a s T s ~ | 2 /3. (20 ) 

1 1 log(V) 
m a = 0, m T6 ~ -^, m 3/2 ~ -, m Ta — . (21) 

with broken SUSY 6 . Here a and r& are the canonically normalized Im(T;,) and Re(Tfe) 
respectively, and ^3/2 = e K / 2 W the gravitino mass; the lightest modulino is the goldstino. 
Thus we can clearly see that the overall volume modulus, r b is the lightest modulus and 
it may have dangerous cosmological effects. The axion a stays ultralight owing to the 
large volume CY, even if e~ 2irTb oc e 2nV is included in the scalar potential. Hence axions 



produced by the modulus decay behave as dark radiation 25|, |26[. By using 01]), the 



modulus partial decay width into the axion pair is calculated as 



1 Tfi 

Tin -»■ aa) = £. (22) 

V ; 48vr M 2 P K J 



5 



Depending on orientifolding and D-brane configuration, moduli is rcnormalized at 1-loop level [3 
Then soft mass on visible branc mentioned below will be changed via such a quantum effect |40(. In 
such cases, we will still have dark radiation via the lightest modulus decay (26 1. 



6 The vacuum has a negative cosmological constant ~ —m^^Mp, however, we will not consider a concrete 
uplifting potential. Including the uplifting potential does not drastically change the results in this 
subsection because of the small negative cosmological constant owing to the SUSY breaking. This 
should be contrasted to the KKLT case. 
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On dimensional grounds, this is an expected order of magnitude, which, therefore, cannot 
be significantly suppressed compared to the other decay processes. Thus there is the 
moduli-induced axion problem, which persists even for heavy modulus mass. In this 
model, the fermionic superpartner of the modulus is the goldstino and obtains a mass of 
m.3/2- Hence the modulus decay into its fermionic component is kinematically forbidden. 
To see how severe the moduli-induced axion problem is, let us consider the modulus 
interaction with the SM sector. 

Suppose that there is a singular cycle moduli denoted by T v which supports the SM 
branes 7 : 

-fr'sin = \l h ^matter (Q, Q ), Wsin = W ma tter(Q), /sin = ~T~ + fo, (23) 

V 47T 

^mat t cr(Q, Q ] ) = e*— /3 [|Q| 2 + (zH u H d + c.c.)]. (24) 

Here Q denotes all visible matter superfields including Higgs ones, fo is a constant and 
we have neglected higher order terms of T v in K mattcv . Vum is the anomalous U{1) gauge 
multiplet, and z is the coupling constant. If Higgs sector has a non-chiral origin like in 



a case of Gauge-Higgs Unification, z ~ 1 is expected [43[; we will consider the case of 
z = 0(1). In addition to (ED), we find 

Tv = 0, m Tv ~ —^. (25) 

Note that T v is absorbed into Vjj(\). In this setup, by using ([7]), we obtain the modulus 
partial decay width to the Higgs boson pair as 

7 Tfl 

T(r b -)• HH) = ^. (26) 

The decay into the higgsino pair is much suppressed due to the approximate no-scale 
structure of the Kahler potential, as mentioned in Sec. [TTJ Thus the branching fraction 
into the axion pair is given by B a ~ 1/(1 + 2z 2 ), which is same as the previous estimate 



7 We can have another odd-parity moduli under the orientifold parity, however it is irrelevant for us 



because it will be also stabilized via another anomalous U(l) D-term, similarly to T v [42 1. 
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(TT4I) because the effective action for Tb possesses a no-scale structure up to a correction 
of 1/V <C 1. Therefore, we need z > 2 or more than four pairs of Higgs doublet to avoid 
the axion overproduction from the modulus decay. 

Let us comment on the soft masses in this model. On the visible branes, the SUSY- 
breaking soft mass is given by 

m 8oft ~ F local ~ ij < m Tb . (27) 

Here, m so ft includes higgsino mass and B-term, and _F local denotes the F-components of 
the local modulus, i.e., dilaton. This is because the SM sector is a local model decoupled 
from the bulk, and the dilaton (and complex structure moduli) are stabilized supersym- 
metrically through a flux compactification, then affected by the a '-correction £/V in the 
Kahler potential; for V = 0(1O 7 ), one finds m soft = 0(1 - 10) TeV and m n = 0(1O 7 ) 
GeV. 8 The decay temperature of the modulus is given by 

If the neutral Wino is the lightest supersymmetric particle (LSP) with the conserved R- 



parity, it will become a good candidate of cold dar 
can produce abundant Winos at m^ ~ 700GeV 



t matter because the modulus decay 



26]. 



B. KKLT: A string-inspired QCD axion model 

Let us again consider a similar effective action on a CY space in flux vacua as a 
generalization of the setup considered in Ref . 15|, ll6| : 



^moduli = -2 log(V); V = (To + Tlf' 2 - ^(71 + T\f 2 - k 2 (T 2 + T 2 f ) 3 / 2 , (29) 

W moduh = W + Ae-^+Be-^ +n ^; a = ^, /3 = ^. (30) 

Here Wq <^C 1 via a fine-tuning of the fluxes, N, M are positive integers, and A ~ B ~ 
0(1). K12 > are constants which depend on intersection numbers among Kahler forms 



We expect that T(, is coupled to the SM gauge bosons through a quantum effect: ^f^rjj-iF^) 2 [44 1. 
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(two cycles) on the CY space, and n £ Z depends on the configuration of an Euclidean 
(instanton) brane. Note that it is only $ = T\ + nT2 that is stabilized in a SUSY way, 
and the axion multiplet A = nT\ — T 2 is absent in the superpotential because we have 
assumed that only the two 4-cycles supported by Re(T ) and Ke(Ti + nT 2 ) are rigid 
divisors. Hence one could not find any other non-perturbative effects 9 . As we will see, 
the imaginary component of A can play a role of the QCD axion to solve the strong CP 
problem. 

At the minimum of the scalar potential, adding a sequestered uplifting potential 5V = 
ee 2if m oduii/3 ^ Q rea ii ze the Minkowski vacuum through a fine-tuning 10 , one obtains solutions 
near the supersymmetric location Dt W ~ Dq>W ~ d^K ~ 0: 

oT^^log^), ReM) . nR e W (z^±|) ; (31 ) 



m„ 



0, m 3/2 ~ W , m s ~ V2m 3/2 , m To ~ m^ ~ log ( — — ) m 3/2 . (32) 



Here, a and s are the canonically normalized Im(^4) and Re(^4.) respectively. We will 
assume that V > is obtained at this minimum, i.e., a < /3. The mixing between the 
saxion and the other moduli is suppressed by a power of l/log(Mp/m 3 / 2 ), hence we will 
ignore such a mixing. Hereafter we will focus on the cosmological effects of s since it is the 
lightest modulus. Note that axino mass is given by m 3 / 2 whereas the modulino masses 
are same as the corresponding moduli masses. Let us estimate the partial decay width of 
saxion s into the axion pair. By using fll]), one finds 



n 3 ^ 2 - n 2 2 Y M 3 S 
7Q8ttk 2 3 Mj 



r» - '- :i::i > a. m 



where 



*-£(^r 5 )>* (34) 



9 For instance, flux can affect the zero mode spectra. However, we will not consider such a possibility 

for simplicity. 
10 For a non-sequestered potential S V non _ seq = ee KlTloduU , we will have similar results [18|, [l9( . 
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Here = Re($). Hence, the saxion partial decay width into the axions is suppressed if 



n^^-f. ,35) 

For n — 1, the axion overproduction is significantly relaxed if the volume of the 4-cycles 
characterized by T\ and T 2 are symmetric under exchange. 

Now let us see the modulus coupling to the SM sector. We will assume that the SM 



localizes on the D-branes wrapping on the cycle supported by T 2 45| : 

1 



-f^matter(V) Q J 



y2/3 



(T 2 + Tl) x \Q\ 2 + (T 2 + Tl) x ^{zH u H d + c.c/ 



(36) 



W vis = H/ matter (Q) , f vis = -^ = — H— ^- . (37) 

An Ait n z + 1 

Here < A < 1. Assuming that the Higgs sector is non-chiral, one will obtain Aqm = A 
and z ~ 1. We will take the minimal case of A = Agm = 1/3 for simplicity. In this 
setup, the imaginary component of A behaves as the QCD axion, since it is massless at 
the perturbative level and obtains a potential dominantly from the QCD instanton effect. 
The decay constant of the QCD axion is given by 



/• = ("' + 1 '^ s W (38) 

As for the soft mass on the visible brane, we will obtain mass spectra through the mirage 
mediation 146. 



F T2 m 3/2 

m soit ~ F ^ 7 ^ < m s . (39) 

T 2 + Tl log (Mr 

6 V m 3, 



If there is the Giudice-Masiero term with z — 1, one finds |/i| 2 ~ Bfi ~ rra 3/2 11 - -^ or 
m 3/2 = £>(10 6 ) GeV, one obtains log(Mp/m 3/2 ) ~ 30 and m soft = C(10)TeV; moduli 
VEVs will be of 0(1 - 10), then f a = C(10 16 ) GeV. 
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In the MSSM, such large (soft) terms in the Higgs sector is dangerous because a color-charge-breaking 
minimum can be realized then 48(. In this paper, however, we will not consider such a problem 
seriously as the Higgs potential depends on the model. 
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The partial decay widths of saxion into the SM sector are estimated from ((Tj) and ffTUl) 
as 

\ 2 K A/f 3 K /If 3 

T(s -> M, WW, ZZ) ~ , 2 ^(2/3)^^, T(s -> A M A M ) ~ AT,^-^; (40) 

On the estimation of the decay fraction into Higgses, we have assumed that the saxion 
decay into the heavy Higgses is not allowed kinematically. The saxion decay rate into 
gauginos is estimated by using 

Bj J* = dAF T> H _L_- , (41) 

which is led from the fact that F Ti / (Tj+ Jf ) is universal for v z. Thus, a sizable contribution 
to that into gauge boson is found from ( fill : 

Note that higgsinos obtain masses comparable to the saxion for z ~ 1, hence decays into 
higgsinos can be forbidden. Otherwise, it would give a comparable contribution to that 
into Higgs bosons. Then the total decay width of the saxion and the decay temperature 
become 

T total c -|- Ml; c ee {n3K ]~f )2 + z 2 sin 2 (2/3)A 2 K 2 + ^, (43) 

487r I0K2 2 

T d K 0.20 GeV x (1 - flj^ (JL) " 2 f «_ ) "' (*) " , (44) 

v ; V30/ \g*{T d )J V 106GeV / 

B a = 2 ^ ~ ^ . (45) 

(n 3 «i 2 - k 2 2 ) 2 + 16/^z 2 sin 2 (2/3)A 2 + Uk\N 9 

There are several ways to suppress the branching ratio into axions. Obviously, if the 
condition ( 135]) is satisfied, the decay rate into axions is suppressed. Alternatively, the 
increase of the decay width into Higgs suppresses B a , which is realized by taking z > 1 
for unsuppressed k 2 . Similar suppression is obtained if there are more than two pairs 
of Higgs doublet. Moreover, for k 2 ^> nK\ the branching fraction becomes suppressed 
due to the relative enhancement of the decay into gauge bosons: B a ~ l/(24N g ). These 
features are seen in Fig. [2] where we plotted the contours of B a on (K2/Ki,n) plane. In 
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FIG. 2: Contours of the saxion branching fraction into the axion pair B a on (k2/ki,ti) plane. 
In this plot we have taken z = and N g = 12. The shaded region is excluded from the axion 
overproduction. 

this plot we have taken z = and N g = 12. The shaded region is excluded from the axion 
overproduction. 

Note that, stnce m „ ft « ™ 3/2 , the LSP S ate likely ovetptodueed by the saxion decay 
into gauginos [37J. To avoid the LSP overproduction, the R-parity should be broken by 
a small amount, and the QCD axion discussed in this section will be a candidate of dark 
matter, if a mild tuning of the initial misalignment angle is allowed. 

For n = 0, T 2 is not stabilized at this level, and hence quantum corrections on the 



Kahler potential and thus SUSY-breaking will stabilize it 20|, |49|, |50(, giving the mass 



only to the saxion Re(T 2 ): 77iRc(T 2 ) ^ m 3/2 whereas mi m (T 2 ) = 0. Thus the saxion will 
mainly decay into axions and the SM gauge bosons; as for the dark radiation the result 
will be similar to the case for n > 0. Whether LSPs can be produced through the decay 
or not depends on the quantum corrections. If the saxion is lighter than the SM gauginos 
it is not necessary for R-parity to be broken, while dark matter mainly consists of QCD 
axion Im(T 2 ) similarly. Otherwise, R-parity should be violated. For n < 0, Kj^ takes 
unphysical values. 
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IV. DISCUSSION AND CONCLUSIONS 

In this paper, we have studied a cosmological fate of the lightest moduli appearing in 
string compactifications, assuming that the moduli dominate over the energy density of 
the Universe. In particular, we have focused on those moduli fields stabilized by SUSY 
breaking effects, as they tend to be lighter than those stabilized in a supersymmetric 
fashion. We have pointed out that those moduli fields can have an important cosmological 
effect, even if they are so heavy that they decay before BBN. This is because the moduli 
often dominantly decay into the light axionic components, which would result in too much 
axion dark radiation in contradiction with observations. Considering that an imaginary 
component of one of the string moduli can play a role of the QCD axion to solve the 
strong CP problem, it is plausible that there is at least one such modulus. Even without 
invoking the QCD axion, the axionic partner of the overall volume modulus is ultra-light 
in the LVS scenario. Such a volume modulus decays into the axion with sizable branching 
fraction. These examples nicely illustrate a new aspect of the cosmological moduli problem 
which cannot be solved just by setting the modulus mass heavy. We have called the axion 
overproduction problem from the modulus decay as the "moduli-induced axion problem" . 

If axions in a string vacuum obtain a small but non-zero mass, the constraint would 
become severer. For instance, they might decay into radiation at late time if sufficiently 
heavy and unstable. They will also contribute to the dark matter and its isocurvature 
perturbations, if stable. 12 To avoid these problems, we may need proper moduli stabiliza- 
tion for the former case, and for the latter case, a fine-tuning of the initial misalignment 



angle or a low-scale inflation will be required additionally [18|, |51 |. 

Now we are in a position to discuss possible solutions to the moduli-induced axion 
problem. The branching fraction can be suppressed in two ways: to increase the partial 



12 If there are multiple axions, one needs to follow their evolution in order of their masses, as the 
isocurvature perturbations depend on their abundance and the size of quantum fluctuations. For 
instance, this will be the case if multiple light axions appear in the SM gauge kinetic function, where 
the QCD axion (mass eigenstate) receives quantum fluctuations of the axions with a possibly different 
decay constant. 



width into the SM sector such as Higgs fields and gauge bosons, or to reduce the partial 
width into the axion pair. For example, the string theoretic QCD saxion can naturally 
have a sizable decay width into gauge bosons. In this case, it is possible that there is a 
small amount of axion dark radiation, which may be detected in the future observations 



52j . On the other hand, the reduction of partial width into the axion pair may imply some 
approximate symmetry on the internal CY geometry. As explicitly studied in Sec. IIHB| if 
there is a symmetry K\ -H- k<i for n — 1, the partial width into the axion pair is suppressed. 
It implies that the size of two 4-cycles are symmetric under exchange. In this case, it is 
also possible that the saxion is stabilized near the low-energy minimum during inflation, 
and its oscillation amplitude can be significantly suppressed. If so, there may be negligible 
amount of axionic dark radiation in the present Universe. Of course, it is always possible 
to assume a huge amount of entropy production by some other (brane) fields which mainly 
decay into the SM particles. 

Depending on models, a heavy field, which is not the superpartner of an axion, also 



can produce the axion dark radiation 
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52J. In this sense, this moduli-induced axion 



problem (and its solutions) lie not only in SUSY models but also non-SUSY models: For 
instance, the inflaton decay in (non-)SUSY models can produce such axion dark radiation, 
if there exist axions in nature. 
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Appendix A: Relevant moduli interactions 

In this appendix we list relevant terms in the Lagrangian for the modulus decay. We 
put several assumptions. First, the modulus T does not appear in the superpotential 
so that its axionic component remains massless. We also assume that T has negligible 
contribution to the SUSY breaking and negligible kinetic mixing between T and the 
SUSY breaking field z: Kt z * = 0. In this Appendix, the canonically normalized fields are 



represented by hats and given by: T = y/K TT T, H ud = JZ ud H ud , A„ = y/Ref^A^. 



X = y/Ref vis X. 



Moduli-axion 



From the axion kinetic term, we easily find 



£ = -±K TTT r(day = -^^rHday. (Al) 



This leads to the modulus partial decay width as (jlj). 

2. Moduli-axino 

From the axino kinetic term we soon find 



£ kin = —K ttt t [~ia)a^ L a + i{d^a ] )a^a] , (A2) 

v2 



which, by using the equation of motion of axino, becomes 

1 lifTT 



--kin 



V2 Kill 



rriaf (aa + h.c), (A3) 



TT 



in terms of the canonically normalized fields. Here the the physical axino mass reads 53] 



ma= K 



1 -^){V T G T + \ Gt G t ) = m 3/2 ( 1 + -^) + %^, (A4) 



where G = K + \nW + InW*, ViGj = d l G j - T^G k , I* = K kl K jTi , and F T = 
—e K ^ 2 K T ^{DjW)* represents the size of modulus F-term. Note that there are in gen- 
eral mixings between the axino and goldstino, but the mixing angle is suppressed by 



20 



~ Kt, which is here assumed to be small. We also assume K tt tKt <C 1. Neglecting 
terms involving Kt, we find the modulus-axino interaction from the axino mass term as 

1 



c 



_ —^ [K TTT (2m 3/2 + F%* + F?*)] raa + h.c. 

2\j2Kr 



(A5) 



, TT 



Noting that Frf* = —m 3 / 2 i n the present approximation, the interaction £ mass vanishes. 
Therefore, the only relevant interaction is £kin, from which we obtain the partial decay 
width as (JSj). 

3. Moduli-Higgs 

The Kahler potential is given by Eq. ($5$). The interaction terms are given by 



\/2 



(d T Z u )r(H u d 2 Hl + h.c.) + {d T Z d )r{H d d 2 H\ + h.c.)] + ^={d 2 r){H u H d + h.c). 

(A6) 



Using the equation of motion for the Higgs boson and modulus, we find 
d T Z u 



£ 



K 



TT 



m H u \Hu\ + ( —^— \ m H d \Hd\ 



9t 

y/2K TT Z u Z c 



m 2 T f(H u H d +h.c). 

(A7) 

Here the Higgs masses are collectively denoted by m 2 Hu and m 2 H . The second term gives 
the partial decay width (|7|). The first terms proportional to OtZ u and dxZ d yields the 
partial width suppressed by the Higgs masses hence we have neglected it. 

4. Moduli-higgsino 



From the higgsino kinetic term, we obtain 



£ kin = V2(d T Z u )r \-iHl<J^H u + t(d^Hl)a»H u 



(A8) 



and similarly for H d . By using the equation of motion of higgsinos, this becomes 

*-kin — 



2 fd T z u 



K 



TT 



m 



/rfH u H d + h.c. 



(A9) 



21 



where /i denotes the physical higgsino mass. Here we have assumed that the higgsino- 
gaugino mixing is negligibly small. 



On the other hand, \x is expressed as 53] 



13 



V- = -7=^e G l 2 V u G d = -±=(gm 3/2 + g T F T "). (A10) 

Then the modulus-higgsino interaction is obtained as 

£ mass = -—=—=== [(2g T + gK T )m 3/2 + g T {F^* + F?*) + 2g TT F T *] fH u H d + h.c. 
V ZK TT Z u Z d 

(All) 
From (1A9I) and (1A11I) . we obtain the partial decay width flHJ). 

5. Moduli-gauge boson 

From the gauge-kinetic function / v j S , we obtain 

£ = -^= (Re(d T f vis )rF^F^ - lm(d T f vis )rF^ a ) 

= ~ A TriAm f V ( R <dTf vis )fF« u F^ - Im(d T f vis )iF« u F^ , 
A^2K TT (Ref vis ) V J 

From this, we obtain the partial width ( TTUj) . 

6. Moduli-gaugino 

From the gaugino kinetic term, we have 

£kin = -£j= [(d T f vis )(i\ a a»d,\ a ) + h.c] , (A13) 

which, by using the equation of motion, is rewritten as 

£ kin = 1 ^Jp-mxf \ a \ a + h.c. (A14) 



^ 



(A12) 



13 If there is an interaction of the Higgs fields with a SUSY breaking field z such as K D z^H u Hd+h..c, it 
would contribute to the [i parameter. However, this does not affect the modulus decay into higgsinos, 
as long as the mixing between T and z is negligible. 
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Here m\ = (dTf V is)F T /(2Ref V i S ) is the physical gaugino mass. From the gaugino mass 
term, we also have 



£ 



{(3r/vta)(*r + F f) + (d 2 T f vis )F T } f\ a \ a + h.c. 



Ay/2K^r(Rel 
From (1A14J) and ( lA15p . we obtain the partial decay width (fTTj) 



. (A15) 
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